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Abstract 

We report on the valence fluctuation of Ce in CeMo2Si2C as studied by means 
of magnetic susceptibility x(T), specific heat C(T), electrical resistivity p(T) and 
x-ray absorption spectroscopy. Powder x-ray diffraction revealed that CeMo2Si2C 
crystallizes in CeCr2Si2C-type layered tetragonal crystal structure (space group 
P4/mmm). The unit cell volume of CeMo2Si2C deviates from the expected lan- 
thanide contraction, indicating non-trivalent state of Ce ions in this compound. The 
observed weak temperature dependence of the magnetic susceptibility and its low 
value indicate that Ce ions are in valence fluctuating state. The formal Lm Ce va- 
lence in CeMo2Si2C <v> = 3.11 as determined from x-ray absorption spectroscopy 
measurement is well bellow the value <z7>~ 3.4 in tetravalent Ce compound Ce02- 
The temperature dependence of specific heat does not show any anomaly down to 
1.8 K which rules out any magnetic ordering in the system. The Sommerfeld coef- 
ficient obtained from the specific heat data is 7 = 23.4 mJ/molK 2 . The electrical 
resistivity follows the T 2 behavior in the low temperature range below 35 K con- 
firming a Fermi liquid behavior. Accordingly both the Kadowaki Wood ratio A/7 2 
and the Sommerfeld Wilson ratio x(0)/7 are in the range expected for Fermi-liquid 
systems. In order to get some information on the electronic states, we calculated 
the band structure within the density functional theory, eventhough this approach 
is not able to treat 4/ electrons accurately. The non-/ electron states crossing the 
Fermi level have mostly Mo Ad character. They provide the states with which the 
4/ sates are strongly hybridized, leading to the intermediate valent state. 
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1 Introduction 



The Ce-based intermetallic compounds have attracted tremendous attention 
with interesting phenomena such as valence fluctuations, heavy fermion be- 
havior, quantum criticality, unconventional superconductivity and Kondo ef- 
fect, etc. [1 ,2.3 .4] These anomalous properties arise due to the hybridization 
between the localized 4/ electrons and conduction electrons. Usually, the com- 
petition between Ruderman-Kittel-Kasuya-Yosida (RKKY) and Kondo inter- 
actions determines the ground state of these compounds. While the RKKY 
interaction favors a long-range magnetic order, the Kondo interactions have 
a tendency to screen the magnetic moments which lead to a non-magnetic 
ground state. [5|6] The result of these two interactions is summarized in the 
Doniach phase diagram. [5] A quantum critical point exists in the regime where 
the strength of RKKY and Kondo interactions are comparable. Quantum crit- 
icality in heavy fermion systems is one of the hot topics in condensed matter 
physics. In Ce-based compounds, quantum criticality and associated unusual 
physical properties such as the non-Fermi liquid behavior are found near the 
magnetic non-magnetic crossover, where the Ce valency is very close to 3+. 
Usually one uses a tuning parameter such as pressure, doping or magnetic 
field to tune the electronic ground state of the system close to quantum crit- 
ical point. Many Ce-based compounds exhibit magnetically mediated uncon- 
ventional superconductivity in the quantum critical regime. Despite the first 
such superconductor, CeCu 2 Si 2 ,[7] was discovered 35 years ago, many fun- 
damental questions have yet not been settled and there is still an intensive 
research going on these superconductors. Pressure experiments on pure and 
Ge-doped CeCu2Si2 revealed two superconducting domes, one at lower pres- 
sure attributed to the magnetic quantum critical point and the other one 
at a higher pressure which could be related to valence fluctuation. [5115] In 
contrast, unconventional superconductivity was not found in Yb-compounds 
near a magnetic quantum critical point, even though there are examples of 
Yb-compounds showing quantum criticality. YbAlB4 is the only compound 
showing superconductivity which is proposed to be connected with critical 
valence fluctuations. [TOlfTTiri2] This provided additional incentive for the in- 
vestigations of valence fluctuation phenomenon. 

In our search for new valence fluctuating compound that may show quantum 
criticality similar to that of YbAlB/t, we have synthesized and investigated the 
physical properties of CeMo 2 Si 2 C . To the best of our knowledge, there is no 
report in the literature on the layered intermetallic compound CeMo 2 Si 2 C. 
The structurally homologous compound, PrMo 2 Si 2 C[i3] has been reported to 
form in CeCr 2 Si 2 C-type ("filled" CeMg 2 Si 2 -type) tetragonal crystal structure 
(space group P^/mmm). The magnetic measurements on PrMo 2 Si 2 C have 
suggested a trivalent state of Pr in this compound. RCr 2 Si 2 C (R = Y, La-Sm, 
Gd-Er), an isoelectronic structural homologue series of CeMo 2 Si 2 C, has been 
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reported to exhibit ferromagnetic ordering of rare earth ions[H] but interme- 
diate valence for R = Ce. |15|[T6] We present the structural and physical prop- 
erties of CeMo2Si2C, viz., x-ray diffraction (XRD), dc magnetic susceptibility, 
electrical resistivity, specific heat and x-ray absorption spectroscopy (XAS) 
and discuss the results in the context of the calculated electronic structure. 
From our measurements, we show that Ce ions are in a valence fluctuating 
state in this compound. 



2 Methods 



2. 1 Experiment 



Polycrystalline sample of CeMo2Si2C has been prepared by the standard arc 
melting technique on a water cooled copper hearth. Since the melting point 
of Molybdenum element is very high, therefore, we first arc melted Si and 
Mo together in the stoichiometric ratio. Then the as-obtained pellet was arc 
melted with carbon and subsequently with cerium. All the constituent ele- 
ments were of high purity (99.9% and above). The samples were flipped after 
each melting and were melted several times to ensure homogeneity. The arc 
melted button was sealed in the Ta crucible and annealed in a vacuum furnace 
at 1500°C for five days. The phase purity of the sample was checked by pow- 
der x-ray diffraction using Cu-K a radiation and scanning electron micrograph 
(SEM). Energy dispersive x-ray (EDX) analysis was used to check the stoi- 
chiometry of the sample. Electrical resistivity measurements were carried out 
using conventional four probe technique in the temperature range 2-300 K in 
a physical property measurement system (PPMS, Quantum design). A com- 
mercial superconducting quantum interference device (SQUID) magnetometer 
(MPMS, Quantum design) was used for magnetic measurements. The specific 
heat was measured by relaxation method in a physical property measurement 
system (PPMS, Quantum design). The absorption spectra of CeMo2Si2C at 
the Ce-L/jj edge (E = 5723 eV) have been recorded in transmission mode at 
the EXAFS beam line Al using the Si (111) double crystal monochromator. 
The absorption behavior has been determined in the energy range 5450 eV 
to 6200 eV with a minimal step size of AE = 0.2 eV close to the absorption 
edge e.g. 5723 ± 25 eV. The powdered sample with particle size smaller than 
20/im was distributed on self adhesive Kapton foil. Multi-layers of 4 foils cor- 
responding to 8 mg/cm 2 resulted an absorption edge step of A/z = 0.5 which 
represents the difference of the low and high energy background functions at 
the absorption edge. Evaluation of spectra, energy calibration and normal- 
ization have been performed with the Athena modul of the Horae software 
package. [17J The simultaneously measured reference spectra of the Ce 3+ com- 
pound CeP04 was used for energy calibration and was measured together with 
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the Ce02 reference material with a dominant Ce 4+ contribution. 
2. 2 Method of calculation 

We have performed the density-functional band structure calculations using 
two full potential codes: local orbital minimum basis band structure scheme 
(FPLO)[[TB] and full-potential linearized augmented plane wave (FLAPW) 
method implemented in WIEN2K code. [12] Perdew-Burke-Ernzerhof (PBE) 
form of the generalized gradient approximation (GGA) was employed for 
the exchange correlation potential. [20J Additionally, to account for the strong 
Coulomb repulsion within the Ce 4/ orbitals a typical value of the Coulomb 
energy U = 6 eV has been chosen for the GGA+U calculations. [2"T1I2"2"] The cal- 
culations were performed using the experimentally obtained lattice parameters 
of CeMo2Si2C where the full-lattice optimizations were done using Rietveld 
refined atomic coordinates. 



3 RESULTS AND DISCUSSION 

3. 1 Experiment 

Figure [1] shows the XRD pattern of CeMo2Si2C together with the structural 
Rietveld refinement profile for the CeCr 2 Si2C-type layered tetragonal crys- 
tal structure (space group PJ^/mmra). Powder x-ray diffraction data were 
collected on the crushed sample and analyzed by Rietveld refinement using 
FullProf software. [23] The obtained lattice parameters for CeMo2Si2C are a 
= 4.213(1) A, and c = 5.376(1) A which are smaller than the values reported 
for homologous PrMo 2 Si 2 C (a = 4.214 A and c = 5.409 A). [13] That means 
CeMo2Si2C unit cell volume does not follow the usual lanthanide contraction 
which is a first indication for a valence fluctuating Ce states. 

The crystallographic parameters obtained from the least square refinement of 
XRD data are listed in Table [H The XRD and SEM reveal the single phase 
nature of the sample. The EDX composition analysis confirmed the desired 
stoichiometry of CeMo 2 Si 2 C. The CeCr 2 Si 2 C-type primitive tetragonal crystal 
structure of CeMo 2 Si 2 C is presented schematically in Fig. [2J This structure 
is different from the very common tetragonal ThCr 2 Si2-type structure: an 
additional carbon atom sharing the same edge with Mo forms a layer of MoC. 
The loss of the body centering results in the c/a ratio 1.27 which is almost 
half of usually observed ratio of about 2.5 in the ThCr 2 Si2-type structure. The 
value of the lattice parameter c in this structure suggests a covalent bonding 
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Fig. 1. (Color online) The powder x-ray diffraction pattern of CeMo2Si2C recorded 
at room temperature. The solid line through the experimental points is the Ri- 
etveld refinement profile calculated for the CeCr2Si2C-type layered tetragonal crys- 
tal structure (space group P4/mmm). The short vertical bars indicate the Bragg 
peak positions. The lowermost curve represents the difference between the experi- 
mental and model results. 




Fig. 2. (Color online) Tetragonal crystal structure of CeMo2Si2C (space group 
P4/mmm). The atomic coordinates are listed in Table I. 

between the Ce-C along [001] direction. [15] The Ce atoms lie at the mid-point 
of the c-axis of the unit cell and there is no direct bonding between them. 

The temperature dependence of magnetic susceptibility x(T) of CeMo 2 Si 2 C 
measured under an applied magnetic field of H = 5 T is shown in Fig. [3] The 
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Table 1 

Crystallographic parameters obtained from the Rietveld refinement of powder x-ray 
diffraction data for CeMo2Si2C. 



Structure 
Space group 
Lattice parameters 
a (A) 

c(A) 

Vcell (A 3 ) 

Refined Atomic Coordinates 



CeCr2Si2C-type Tetragonal 
P4/mmm 

4.213(1) 
5.376(1) 
95.42(1) 



Atom 


Wyckoff 
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Ce 


lb 


0.0000 


0.0000 


0.5000 


Mo 


2/ 


0.0000 


0.5000 


0.0000 


Si 


2h 


0.5000 


0.5000 


0.2812 


C 


la 


0.0000 


0.0000 


0.0000 



very low value of susceptibility and its weak temperature dependence down 
to 100 K clearly indicate that Ce ions are in valence fluctuating state. Similar 
shape of susceptibility curves has been observed and are typical for Ce-based 
valence fluctuating compounds. [24T251I26II27] Sometimes the magnetic suscep- 
tibility does show a curie tail at lower temperatures in valence fluctuating sys- 
tems for instance in CeoCogGe^.p^ Ce 2 Ni 3 Si5,[27] and CeNi 2 B 2 C. (28] which 
arises due to free Ce 3+ magnetic impurity ions. 

The dc magnetic susceptibility for a valence fluctuating system can be de- 
scribed with the two level ionic interconfiguration fluctuation (ICF) model. 
The theory of ICF was first proposed by Hirst [30] and latter Sales and 
Wohlleben developed it further to explain the valence fluctuating behavior 
observed in few Yb-based compounds. [31] According to the ICF model. [32] 
the temperature dependence of the susceptibility is given by 
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Fig. 3. (Color online) Temperature dependence of the magnetic susceptibility of 
CeMo2Si2C at 5 T. The solid red line through the magnetic susceptibility data 
shows a fit to the two-level ionic interconfiguration fluctuations (ICF) model in the 
temperature range 10-300 K. 

and 



T 2 + T 2 sf 



1/2 



(3) 



where /i n and the effective moments in Af n and 4/( n ^ states, and 

(2J n + 1) and (2J n _x + 1) are the degeneracies of the corresponding energy 
states E n and E n _i. Here E ex is interconfigurational excitation energy which 
is equals to (E n -E„,_i) and T s j is the spin fluctuation temperature associated 
with the valence fluctuation. In order to take care of the contribution of Ce 3+ 
ions which belong to the impurity phases present in CeMo2Si2C, we have added 
the term Xim P = n * C/(T — 9) to the magnetic susceptibility. Thus the final 
equation is 
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with 



Z/(T) 1 + 6exp{-E ex /K B T*) (5) 

where Xo is temperature independent term and n is the fraction of stable Ce 3+ 
ions. Here we have taken Ce 4+ (J = and \i = fig) state as ground state 
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and Ce 3+ (J = | and /i = 2.54 /x B ) as an excited state. The solid line through 
the magnetic susceptibility data in Fig. [3] is the fit to the equation 4. The 
values of parameters obtained from the least square fits are n = 0.06, E ex = 
677 ± 12 K, T aJ = 205 ± 3 K, 9 = 1.9 K and xo = -1-67 x 10" 4 emu/mole. 
The values of E ex and T s f are in the range of typical values found for the 
valence-fluctuating system Ce 2 Rh 3 Si 5 (E ex = 845 K and T s f = 129 K).[22] 
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Fig. 4. (Color online) Temperature dependence of the magnetic susceptibility of 
CeMo2Si2C at 5 T. The solid line through the magnetic susceptibility data in the 
temperature range 15-80 K shows a fit to the equation x(T) = xo+ n *C/ (T—9). The 
blue circles in the inset represent the susceptibility data obtained after subtracting 
the extrapolated impurity component from the measured susceptibility data in the 
temperaure range 10-300 K and the solid line through the data points represents 
an ICF fit as discussed in the text. 



Furthermore, to understand the intrinsic behavior of the susceptibility for 
the sample i.e. xCO without the contribution from Ce 3+ impurity phase, 
we have subtracted the magnetic impurity contribution from the measured 
susceptibility data. The susceptibility data were first fitted to the equation 
x(T) = Xo + n * C/(T — 9) in the temperature range 15-80 K where the domi- 
nant magnetic contribution is coming from Ce 3+ magnetic impurity pase (see 
Fig. 4). The parameters obtained from the least square fit, n = 0.06 and 9 = 
2.32 K are nearly same as found in the impurity part of ICF fitting discussed 
earlier. In next step, the impurity part Xim P = n* C/(T — 9) was subtracted 
from the measured x{T) in the temperature range 10-300 K. The resulted 
(after subtraction) susceptibility x ~ Ximp posses a broad maximum which is a 
characteristic feature of a valence fluctuating system [221133] (see inset of Fig. 
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4). The 

X ~ Ximp data were fitted to the following ICF equation 
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where the impurity concentration n is kept fixed at 0.06. The fit to the data 
yields the values E ex = 701 ± 4 K, T sf = 217 ± 1 K and xo = -1-24 x 
10~ 4 emu/mole, which are close to the values obtained from the ICF fitting 
earlier above. 




Fig. 5. (Color online) The temperature dependence of zero field-cooled (ZFC) and 
field-cooled (FC) dc magnetic susceptibility in the temperature range 2-40 K mea- 
sured under various applied magnetic fields. 



Even though the valence fluctuating systems are not expected to order magnet- 
ically, the magnetic susceptibility measured under low magnetic fields exhibit 
sharp increase below 8 K (Fig. [5]). However, the low value of magnetization 
at 5 K (~ 0.003 /is/Ce at 0.01 T) and the absence of any anomaly in the 
specific heat data suggest that this low temperature magnetic susceptibility 
anomaly is not intrinsic and arises from small amount of impurity phase. The 
field dependence indicate this impurity phase to be ferromagnetic. Therefore 
it is likely a small amount of CeSi2-x-[34j 

The specific heat C(T) of CeMo 2 Si 2 C measured under zero applied magnetic 
field is presented in Fig. [6] in the temperature range 1.8-35 K. The specific 
heat exhibits no anomaly in this temperature range down to 1.8 K and hence 
confirms the absence of magnetic ordering in CeMo2Si2C. We could fit the low 
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Fig. 6. (Color online) Temperature dependence of the specific heat data, C(T) 
of CeMo2Si2C measured in zero field in the temperature range (1.8-35 K). The 
solid line through the specific heat data in the inset shows the fit to the equation 
C(T) = 7 T + /3T 3 . 

temperature specific heat data below 12 K to the equation 

C(T) = 7 T + (3T 3 (7) 

where 7T is the electronic contribution to the specific heat, /3T 3 is the phononic 
contribution to the specific heat (inset of the Fig. [6]). The fit to the data 
yields the value of 7 = 23.4 mJ/mol K 2 and (3 = 1.6 x 1CT 4 J/mol K 4 . 
This gamma value is at the lower border of the range expected for valence 
fluctuating Ce system, but comparable to the value found e.g. in CeRu 2 (7 = 
23 mJ/mol K 2 ),[35j Ce 2 Co 3 Ge 5 (7 = 17 mJ/mol K 2 )[24j and CeRu 3 Si 2 (7 = 
39 mJ/mol K 2 )[29]. The Debye temperature Qd for CeMo2Si2C was estimated 
from (3 using the relation [36] 

/l27r 4 niA V3 

HH (8) 



where R is molar gas constant and n = 6 is the number of atoms per for- 
mula unit (f.u.). Therefore, we obtain Debye temperature 0£> = 417 K for 
CeMo 2 Si 2 C. 

The temperature dependence of electrical resistivity p(T) of CeMo 2 Si 2 C is 
shown in Fig. [71 The electrical resistivity data exhibit metallic behavior with 
residual resistivity, p ~ 12.7 fifl cm at 2 K. The large value of residual 
resistivity ratio (RRR = P3ook/P2k) ~ 9 indicates good crystallinity of our 
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Fig. 7. (Color online) Temperature dependence of electrical resistivity data, p(T) of 
CeMo2Si2C as a function of temperature. The solid line in the inset shows the fit 
to the equation p(T) = po + AT 2 . 

sample. We do not observe any prominent anomaly in the resistivity data 
and the data fit well to the following equation in the low temperature range 



The fitting yield the parameters po~12.3 fi f2 cm for the residual resistivity 
and the coefficient A = 2.57 x 10~ 3 p Q cm K -2 . Thus, the electrical resis- 
tivity of CeMo2Si2C exhibits the Fermi liquid behavior as usually observed 
for Ce-based valence fluctuating systems. |37] Within the Fermi liquid theory 
the A coefficient is related to r y 2 . Kadowaki and Woods found that in heavy 
fermion systems and valence fluctuating systems the ratio A/7 2 is of the order 
of 1 x 10~ 5 p Q cm (mol K/mJ) 2 .[38j Later on N. Tsujii et al.[3H] suggested 
this ratio to scale with 2/N(N-l) where N is the degeneracy of the orbital state 
of the / element. For an intermediate valent Ce system N = 6 gives A/7 2 = 
6.7 x 10~ 7 fi Q cm (mol K/mJ) 2 . For CeMo2Si2C our experimental data result 
in A/7 2 = 0.5 x 10~ 5 fi Q cm (mol K/mJ) 2 , just the original Kadowaki Wood 
value. Furthermore, in a Fermi liquid one expect a T-independent susceptibil- 
ity Xfl at low temperatures with a value scaling with 7 too. This is expressed 
in the Wilson Sommerfeld ratio i.e. 



which is expected to be close to one. Taking the low temperature value of 
Xfl~ 5.6 x 10~ 4 emu/mole from the impurity corrected x(T) shown in the 
inset of Fig.4, 7 = 23.4 mJ/mol K 2 as determined above, and fi e ff = 2.54 /ig as 
expected for the free Ce 3+ moment we get Rw = 0.81, while taking p e ff = 1.73 



10-35 K , i.e. 



p{T) = po + AT 2 
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as for a free conduction electron, one get Rw = 1.7. Thus, low T susceptibility, 
specific heat and resistivity data fulfill Fermi liquid expectations. 
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Fig. 8. (Color online) L///-edge x-ray absorption of CeMo2Si 2 C at 300 K. The solid 
black line through data points shows the least square fit. The contributions due to 
4/ 1 and 4/° configurations are represented by the solid Gaussian lines, while the 
dashed line represents the arc tan function. 

The Ce-Luj-edge x-ray absorption spectrum of CeMo2Si2C taken at room 
temperature is presented in Fig. [HJ To analyze the spectrum in detail two 
Gaussian functions were taken into consideration which represents the 4/ 1 
and 4/° configuration. Additionally, an arc tan function was added to provide 
the relative weight of the two electron configurations. The larger peak at 
5724 eV corresponds to the 4/ 1 configuration and the smaller peak at 5734 eV 
corresponds to the 4/° configuration. The Ce-Lm clearly indicates that the Ce 
ions in this compound are in valence fluctuating state. Moreover, the average 
formal Ljjj valence of Ce can be obtained from the intensity ratio of the two 
Gaussian peaks using the equation <u> = 3+Ii/(Ii+I 2 ). Hence we get the 
value of average formal Lni valence of Ce atoms in CeMo2Si2C as <V> = 3.11 
at room temperature. This value is less than the value (~ 3.4) obtained for 
a tetravalent CeO2 440j Similar values of mean formal Lni valence of Ce ions 
have also been reported for other valence fluctuating compounds, for example, 
CeRhSi 2 (<V> = 3.15), Ce 2 Rh 3 Si 5 (<V> = 3.18)[22j and CeNi 2 B 2 C (<u> = 
3.18). [28J at room temperature. 



3.2 Band Structure Calculations 



In order to gain some information on the electronic state, electronic structure 
calculations were performed for the compound using two different approaches, 



12 



> 



CO 
O 
Q 



40 
30 
20 
10- 





CeMo 2 Si 2 C ; 






total DOS GGA \ 






total DOS GGA+U 1 

1 




Ce 5p 

1 

f 


I 
■ 
1 

i 1 i 1 i 1 i 1 i 1 i ■ 





2.0 
1.5 
1.0 
0.5 
0.0 
1.5 
^ 1-0 

e" 

o 

■£ 0.5 
> 

CD 

"55 0.0 

£ 

"S 

W 1.5 
to 

8-i.o 

0.5 
0.0 
1.5 
1.0 
0.5 
0.0 



-15 -12 -9 -6 -3 
Binding energy (eV) 



GGA+U 

- . Ce 5p 
...Ce 4d 
— Ce 4f 


- ■i < ^*~^ V | 






- • Mo 5s 
• ■ - Mo4p 

— Mo 4d * 






--- Si 3s : 
— Si 3p : 

i > 
■■ ■ 

.A. , ^/Wj 




1 1 — 1 

...C2s 
— C 2p 


"1 [ 





-12 



-6 -3 
Binding Energy (eV) 



Fig. 9. (Color online) Total (upper panel) and partial densities of states (bottom 
panels) for CeMo2Si2C. Here the Fermi level (^f) corresponds to zero binding 
energy. 

i.e., with FPLO and full potential LAPW method. Both the methods give very 
similar results and therefore only those obtained by using the FLAPW method 
are discussed here. Because of the strong correlations in the 4/ shell, the den- 
sity functional approach is not able to treat the 4/-electrons of Ce accurately. 
In contrast non-/ valence states are usually only weakly affected by these cor- 
relations, except close to the Fermi level where hybridization with /-states 
becomes relevant. Therefore, such a calculation allows insight into the non-/ 
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valence states and gives some information on the states which are hybridized 
with the Ce 4/ states at the Fermi level. In this method the unit cell is divided 
into non-overlapping muffin-tin spheres centered at the atomic sites and an 
interstitial region. The muffin-tin radii were chosen to be of 2.5 a.u. for Ce, 
2.1 a.u. for Mo, 2.2 a.u. for Si and 1.86 a.u for C. The set of plane-wave expan- 
sion Km ax was determined as Rmt^Kmax equal to 7 and K mesh used was 
13 x 13 x 10. We have performed the calculations with and without spin polar- 
ization. Notably the self-consistent results for magnetic solutions were nearly 
the same as those for non-magnetic ones because the spin polarized calcula- 
tion resulted in a non-magnetic state. Total and local magnetic moments were 
approximately zero (~ 10~ 5 /iB/Ce) in the spin polarized calculations. 

Figure [9] shows the total and partial Density of States (DOSs) for CeMo 2 Si 2 C, 
calculated within the GGA and GGA+U (U = 6 eV) approximations. To 
describe in detail, here we divide the total DOS into four main parts. 

(1) The peak located below -16 eV is due to Ce 5p electronic state. It is re- 
markable that the GGA+U peak is shifted towards higher binding energy 
by 1 eV with respect to GGA. 

(2) The sub-band situated between -12 eV and -6 eV, is formed mainly by 
Si 3s and C 2s states. 

(3) The third part, which is located between the -6 to -2 eV contains strongly 
hybridized Mo Ad , Si 3p and C 2p states. Hybridized states of Ce 4d and 
C 2p are observed as well just below the Fermi level (see the partial 
density of states of Fig. 

(4) The bottom of the conduction band or at the Fermi level, the main con- 
tribution to the DOS is provided by Ce 4/ and Mo Ad electronic states. 
Hybridization between Ce 4/ and Mo Ad states are observed near the 
Fermi level. In case of GGA+U the main peak (Ce 4/) is broadened and 
slightly shifted towards right side above the Fermi level in comparison to 
GGA and consequently reduced the occupancy of Ce 4/ electrons at the 
Fermi level. 

It should be noted that at the Fermi level a small contribution to the total 
DOS come from Ce 5d, Si 3p and C 2p states as well. The orbital occupancy of 
electrons inside the atomic sphere of the Ce atom is calculated by the GGA ap- 
proach as 6s L9S 5d a64 4/°- 94 and by GGA+U calculation as 6s L " 5rf a85 4/ a64 . 
Thus the f-count in these calculations is slightly (GGA) or far (GGA+U) be- 
low 1. Comparison with the LI 1 1 result shows that GGA overestimate the / 
count while GGA+U provides more realistic value of / electron count. 

Based on the DOS at the Fermi level here we obtained the Sommerfeld coeffi- 
cients 70 and the Pauli paramagnetic susceptibility xfl as calculated under the 
assumption of free electron model, 70 = (vr 2 /3)N[Ef\K^ and \fl = HbN[Ef], 
where N[Ep] is the electron density of states at the Fermi level. All the calcu- 
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Table 2 

Total densities of states N[Ep] (in states/eV/f.u) and partial densities of states 
(in states/eV/atom) at the Fermi level, the calculated Sommerfeld coefficient 70 (in 
mJ/mol K 2 ) and molar Pauli paramagnetic susceptibility xfl (in 10 emu/mol). 



Method 


Total N[E F ] 


7o 


X 


Ce-4/ 


Mo-4d 


GGA 


3.86 


9.10 


1.24 


1.47 


0.52 


GGA+U 


3.16 


7.45 


1.01 


0.77 


0.52 



lated values are listed in Table |2j The calculated 70 and xfl values are a factor 
of about 3 and 6, respectively, below the experimental values. This difference 
can be very simply attributed to the strong renormalization connected with 
the valence fluctuations. 




Fig. 10. (Color online) Band structure of CeMo2Si2C along the symmetry lines of 
the Primitive tetragonal BZ. Mo Ad contribution to the bands is represented by the 
larger circles. Different colors of the bands are used just to guide the eyes. 

The band structure along various symmetry directions calculated within the 
GGA is shown in Fig. [KB The band structure near the Fermi level is mainly 
composed of three bands which are crossing Ep- The fat part (bigger circles) of 
the bands represent the contribution from Mo 4d states which are considerably 
hybridized with the Ce 4/ bands near E^. Two hole like bands cut by the Fermi 
level near T and X points and a highly dispersed electron like band crosses the 
Fermi level near M and A points. The strong dispersion along T-Z shows that 
this compound, despite having a quasi 2D structure, has a true 3D electronic 
structure. Interesting features are observed between T-Z and A-Z as the sharp 
tip of the bottom of the conduction band touches the sharp tip of the top of 
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the valence band. 



4 Conclusion 

To summarize, our comprehensive studies of physical properties of CeMo2Si2C 
by x(T), C(T), p(T), and XAS measurements reveal the valence fluctuating 
state of Ce in CeMo2Si2C. The compound crystallizes in the CeCr 2 Si2C-type 
layered tetragonal structure (space group P^/mmm). The value of magnetic 
susceptibility is very low and almost independent of temperature except a 
curie tail at lower temperature (as observed in many valence fluctuating com- 
pounds) which arises due to a small amount (~ 6 %) of an impurity Ce 3+ 
phase present in the sample. An ICF analysis of x(T) result in a fluctuation 
temperature T s f = 205 ± 3 K and an interconfigurational excitation energy 
E ex = 677 ± 12 K. The temperature dependence of specific heat at zero mag- 
netic field does not show any anomaly throughout the temperature range and 
hence confirms the absence of any magnetic ordering down to 2 K. The Som- 
merfeld coefficient obtained from the specific heat data 7 = 23.4 mJ/mol K 2 
is relatively low for a valence fluctuating compound which indicates a strong 
hybridization of Ce Af states with the conduction states at the Fermi level. 
The electrical resistivity exhibits a metallic behavior with a T 2 dependence in 
the low temperature range indicating the Fermi liquid behavior as observed for 
many Ce-based valence fluctuating systems. XAS studies indicate an average 
formal L IH valence <V> = 3.11 hence gives a direct evidence for valence fluc- 
tuations. Nevertheless, this value is less than that obtained for a tetravalent 
CeC>2 (<v> ~ 3.4). The density functional calculations show strong disper- 
sions along all direction evidencing a 3D electronic structure despite the quasi 
2D chemical structure. Strong hybridization is found between Mo Ad, Si 3p 
and C 2p states, and between Ce 4d and C 2p states at lower binding ener- 
gies below the Fermi level. The bands at the Fermi level bear a large Mo Ad 
character and provide the states which are hybridized with the Ce Af states. 
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